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We find that multifractal scaling is a robust property of a large class of continuous stochastic 
processes, constructed as exponentials of long-memory processes. The long memory is characterized 
by a power law kernel with tail exponent ip + 1/2, where ip > 0. This generalizes previous studies 
performed only with tp = (with a truncation at an integral scale), by showing that multifrac- 
tality holds over a remarkably large range of dimensionless scales for <p > 0. The intermittency 
multifractal coefficient can be tuned continuously as a function of the deviation ip from 1 /2 and of 
another parameter a 2 embodying information on the short-range amplitude of the memory kernel, 
the ultra-violet cut-off ("viscous") scale and the variance of the white-noise innovations. In these 
processes, both a viscous scale and an integral scale naturally appear, bracketing the "inertial" scal- 
ing regime. We exhibit a surprisingly good collapse of the multifractal spectra £(g) on a universal 
scaling function, which enables us to derive high-order multifractal exponents from the small-order 
values and also obtain a given multifractal spectrum £(q) by different combinations of (p and a 2 . 

PACS numbers: 



I. INTRODUCTION 



Generalizing the cascade models that started with 
Richardson [lj and Kolmogorov 2], multifractal cas- 
cades have been introduced in turbulence 0,01 to model 
the anomalous scaling exhibited by the moments of the 
velocity increments in hydrodynamic turbulence (see for 
instance p| and references therein). They have been 
since applied to many other complex fields including 
fractal growth processes, geophysical fields, high energy 
particle physics, astronomy, biology and finance 6]. The 
constructions involved in multifractal cascades are based 
on hierarchical geometries coupled with multiplicative 
noise and form discrete hierarchical cascades They 
have been very useful to highlight a general mechanism 
for intermittency and multifractality which reflects the 
presence of intermittent bursts of fluctuations with long- 
range correlations. Accordingly, the long-range correla- 
tions are seen to result from the large-scale structures 
that impact the smaller scales through a hierarchical 
cascade. But discrete cascades have limitations and de- 
fects such as spurious effects due to the discreteness 
(scaling holds only for certain scale ratios 8J), non- 
stationarity, and absence of causality in the time domain 
(see however 
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Here, we study the multifractal properties of a class 
of continuous stochastic processes, constructed as expo- 
nentials of long-memory processes with power law mem- 
ory. Previous works briefly reviewed below have been 
concerned with the case where the power law memory 
has a tail exponent equal to 1/2, which leads to a log- 
arithmically decaying correlation functions and the ne- 
cessity for a regularization at large time scales, i.e., the 
introduction of a so-called integral scale. With logarith- 
mic correlation functions, previous works have shown 
that this class of processes exhibit the property of mul- 
tifractality. Here, we extend the problem to power law 
memory with tail exponent ip + 1/2 which can vary ar- 
bitrarily above 1/2. As a consequence of the faster- 
than-logarithmic decay of the correlation function of 
the process, the property of multifractality cannot hold 
exactly anymore. We show however that multifractal- 
ity holds over a remarkably large range of dimension- 
less scales and that the intermittency coefficients can 
be tuned continuously as a function of the deviation ip 
from 1/2 of the exponent of the power law memory and 
of another parameter a 2 embodying information on the 
short-range amplitude of the memory kernel, the ultra- 
violet cut-off scale and the variance of the white-noise 
innovations. For this, we present a motivated robust 
algorithm to determine the exponents of the mul- 
tifractal spectrum that we apply on our numerically de- 
termined structure functions or moments of the stochas- 
tic process. We exhibit a surprisingly good collapse of 
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the multifractal spectra on a universal scaling function, 
which enables us to derive high-order multifractal expo- 
nents from the small-order values. The scaling ansatz is 
validated by direct numerical evaluations of integral ex- 
pressions of the moments of orders up to q — 5. Our 
results offer an interesting generalization of the class 
of multifractal random walks introduced recently and 
provide a physically interpretable source of multifrac- 
tal intermittency in terms of the parameters cp and a 2 . 
Our results have potential use in all the fields in which 
multifractal properties have been discussed in the time 
domain. For instance, they provide a rational for the 
approximate multifractal signatures observed in simple 
agent-based models of social networks without the 
need to justify an exact logarithmic scaling for the cor- 
relation function of the logarithm of the observable. 

The organization of the paper is the following. Sec- 
tion 2 presents a short review to position the stochastic 
process which we also define. Section 3 defines the ef- 
fective multifractal exponents and present a simple ef- 
ficient method for their estimation. Section 4 focuses 
on the scaling properties of the second-order moment 
and its associated exponent £(2). Section 5 presents our 
main results for the higher-order moments and their ex- 
ponents ((q). Section 6 offers some numerical and math- 
ematical insights on the origin of the effective multifrac- 
tality. Section 7 concludes. 



II. STOCHASTIC CONTINUOUS PROCESSES 
AS EXPONENTIALS OF PROCESSES WITH 
POWER LAW MEMORY WITH ARBITRARY 
EXPONENTS 

Recently, inspired by the logic of the construction 
of discrete hierarchical cascades, several works have di- 
vised genuine stochastic continuous stationary processes 
which reproduce their main properties 0, E2, EH EH 
EU EH E3 ■ In particular, the so-called multifractal ran- 
dom walk (MRW) has been introduced by Bacry, Delour 
and Muzy as the only continuous stochastic stationary 
causal process with exact multifractal properties and 
Gaussian infinitesimal increments Sornette et al. 

|15| have shown that the increments 5 T X(t) at finite 
scale t of the MRW can be approximated by 



S T X(t)= f 



e(t) 



r(t) 



where 8 T X(t) = X(t + r) — X(t), e(t) is a standardized 
Gaussian white noise independent of u> T (t) and lo t (t) can 
be expressed as an auto-regressive process 



Jr (t) = 



dW(t') h T (t-t') 



(2) 



where W{t) denotes a Wiener process with unit diffu- 
sion coefficient and the memory kernel h T (-) is a causal 
function specified by its Fourier transform 



[h T {f)f = 2A 2 r 1 



Tf 



sin(i) 



dt + 0(fT]n(fr)) 



The expression @ shows that 



(3) 



h T {t') ~ K 



X 2 T 



for r < t' < T , (4) 



where the so-called integral scale T delineates the 
boundary beyond which the correlation vanishes exactly. 
This slow inverse square-root power law decay (@J of 
the memory kernel in J2J) ensures the long-range loga- 
rithmic dependence of the correlation function of lu t (t) 
p"H , which is one important ing redient for the multi- 
fractality of 5 T X(t). Schmitt has studied in details 
the stochastic process 111 12ft with a kernel h T (t') exactly 
given by the square-root power law (0J for t' between 
scale t — 1 and T, with a smooth regularization for 
t < t. Not surprisingly, this process exhibits multifrac- 
tality in the range of scales between r and T. 

We study the positive stochastic process S T X(t) gen- 
eralizing 1)1121) with 0}, defined by 



S T X(t) 



^{t')dt' , with fj,(t) = k e" (t) , (5) 



with ui(t) of the form 10) 



dW(t') h(t - t') , 



(G) 



and 



h(t) 



h 



(1 + 2^+1/2 



H(t) , x = t/t, (7) 



where t is some "microscopic" characteristic scale, reg- 
ularizing the singularity of the power law in the propa- 
gator h(t) at t = and H(t) is the unit step Heaviside 
function ensuring the condition of causality inherent to 
most applications. The main departure from the pre- 
viously cited works is to consider an exponent h + <p 
larger that 1/2 for the power law decay of the memory 
kernel h(t). As already mentioned, one can prove rigor- 
ously [Tj| that the process (|f!2[) with 0, whose power 
law approximation for the memory kernel has the exact 
exponent 1/2, i.e. ip = 0, has a logarithmic decaying 
covariance function of the auxiliary stationary Gaussian 
stochastic process w(t) 



(uj(t)uj(t + t)) 



In [ - 



T < T 



(8) 
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associated with the multifractal signature of the process 
X(t) given by 

([6 T X(t)] q ) = a(q)r c ^ , for r < T . (9) 

In expression © , the angle brackets denote a statistical 
averaging and the multifractal "spectrum" has the 
parabolic form 



(10) 



where A 2 = — C"(0) is the so-called intermittency co- 
efficient. In contrast, the existence of multifractality 
defined by the nonlinear spectrum £(q) has not been 
studied previously for the process (|5I7|I with non-zero 
values ip > 0. It is not obvious a priori that multifrac- 
tality will be observed because the deviation from 1/2 
for the exponent of the power law decay of the mem- 
ory kernel implies that the covariance function of the 
stochastic process u>(t) is no more logarithmic, which 
was the fundamental reason for the existence of multi- 
fractality in the MRW. However, Ouillon and Sornette 
[l8L Il9l | have recently shown that the process 1)51711 with 
non-zero values tp > has robust multi-scaling proper- 
ties. They have derived this process from the physics of 
thermally activated rupture and long memory stress re- 
laxation for earthquakes, and have shown that this pro- 
cess predicts that seismic decay rates after mainshocks 
follow approximately the Omori law ~ l/t p with expo- 
nents p(M) linearly increasing with the magnitude M of 
the mainshock, in agreement with observations pHIT^. 
Such multi-scaling suggest that the property of multi- 
fractality in the sense of @ should also be present. 

A significant difference between the process (|5I7|) with 
non-zero values tp > and with tp = is that no integral 
scale T is needed to regularize the theory. Furthermore, 
all moments ([S T X(t)] q ) are finite for tp > whereas all 
the moments of order q > q* for ip = 0, where g» satis- 
fies to equation C(<?*) = 1, arc infinite. This divergence 
signals that the probability density function (PDF) of 
the increments S T X(t) is heavy-tailed with an exponent 
smaller than q* [2(j. The absence of divergence of all 
moments of the PDF for tp > excludes a heavy tail but 
not fat tails of the PDF of increments, such as stretched 
exponentials (which are known to approach arbitrary 
well any power law, see chapter 2 of 21]). Thus, if mul- 
tifractality exists for ip > 0, it may be observed for any 
q > 0. 

In addition to tp, the other key parameter controlling 
the multifractal properties of the process (1~>I7I will be 
shown to be 



a 2 = 



h 2 (t)dt = hl — 



(11) 



integral scale must be re-introduced to regularize the 
theory. 



III. DEFINITION AND DETERMINATION OF 
EFFECTIVE MULTIFRACTAL EXPONENTS 

A. Definitions and notations 

A general theoretical characterization of the incre- 
ments (JjjJ is offered by the moment functions defined 
by 

i q 
M(t 1 ,...,t q ) = (]J^t r )) = (^) J] G(x k -Xi) 

k=i + l 

(12) 

of the lognormal density fj,(t) = ft defined in JSJ. 
In 112|) . we use the following notations: Xi — U/i, 



9\ = 



G(y) = e- CTd ^, d(y) = l-C(y), 

(13) 

where C(t) is the normalized (C(0) = 1) covariance 
function of the Gaussian process uj(t) defined in JBJ, 



C 



(j) =^ J o °°h(t)h(t + r)dt. (14) 



We start our investigation of the multifractal proper- 
ties of the increments S T X(t) defined in (jjjj by calculat- 
ing the moment (\8 T X{t)] q ) for q — 2 and checking if 
the power law scaling of the form holds with an ex- 
ponent C(2) smaller than 2. If this is the case, and from 
the fact that £(0) = and £(1) = 1, we can conclude 
that C(g) is a nonlinear function of q, the hallmark of 
multifractality. The fact that £(1) = 1 results from the 
positivity of the stationary density (i(t). 

In order to get the exponent C(2), let us study the 
normalized second moment 



= pT^^rX{t)f) = J\ Xl J\x 2 G{x 2 -x x ) . 

(15) 

For the numerical analysis of l|15fl . we use the more con- 
venient representation 



S 2 {y) = 2 / (y-x)G(x)dx 
Jo 



(16) 



The question we address is whether and how 5*2(2/) can 
be approximated by the power law 



Note the divergence of a 2 for tp — > 0, for which the 



S2(y) = A 2 y^ . 



(17) 
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B. Properties associated with the second-order 
moment 

First, let us determine the conditions on G(x) for 
which the form (|17f) is exact. For arbitrary behaviors 
°f ^2(2/)) one can always introduce a local exponent de- 
fined by 



C(2,y) = 



dlny 



(18) 



which recovers C,(2, y) = £(2) if the power law l|17fl holds 
exactly. In the general case where £2(2/) is no t a pure 
power law, we can write 



C(2,y) 



A(y) 



where 



A(y) 



2X(y) 



xG(x)dx 

y 

G(x)dx 



(19) 



(20) 



Thus, £(2, j/) is independent of y and the scaling law (|17fl 
is exact if X(y) is proportional to y, where X{y) has the 
interesting interpretation of being the barycenter of the 
segment [0, y] whose mass density is G(y). 

It is easy to show that X(y) ~ y if and only if 
G(y) is a constant or a pure power law G(y) ~ 2/ A , 
with some exponent that we denote < A 2 < 1 for 
a reason that will be clear soon. In the former case, 
X{y) — y/2, A(y) = 1, which yields the non- fractal 
scaling ((2,y) — 2. From fT^I . we see that G(y) is a 
constant if the covariance function C(t) of the Gaus- 
sian process u>(t) defined in (|14|) is also a constant that 
is dC(r)/dT = 0. Taking the derivative of i|14|) with re- 
spect to r, integrating by part and equating to zero for 
arbitrary values of r imposes the condition h(0) = 0, 
which just means that the measure n(t) is uniform, 
hence the trivial exponent £(2, y) — 2. The other case 
G(y) ~ y~ x is more interesting. This yields 



X(y) 



1 



2- A 



A 2 

— V 



1-A 2 



(21) 

From (|13[) . we see that this case corresponds to C(r) 
being an exact logarithmic function of t, which is the 
property already mentioned above with JSJl at the origin 
of the exact multifractality of the process 6 T X(t). 

This discussion shows that, for <p > 0, the scaling 
(|17(l of the second-order moment can only hold approxi- 
mately at best. In particular, the local exponent £(2, y) 
has a simple regular behavior at the two boundaries 
y — > and y — > 00. Indeed, due to the limits 



lim G{x) = 1 



lim G(z) - (/i) 2 /</i 2 



(22) 



we have 

S 2 (y)^2/ 2 (1,-0), 5 2 (y)^^y 2 (y-oo), 



so that 



</< 2 



hm C(2,y)= lim C(2,y)=2 



(23) 



(24) 



Borrowing the terminology from hydrodynamic turbu- 
lence, these two limits ()24|l implies the existence of an 
effective "viscous" scale t v and of an integral scale Tj, 
such that if r < r„ and r > Tj, then the increments 
JSJ are not multifractal. Note that the limit J/ — 00 is 
attained by taking I — at fixed t since y = t/l. The 
limit I — thus recovers a trivial regime, in contrast 
with the MRW in which the role of £ is played by the 
scale of resolution which, when going to zero yields a 
non-trivial genuine multifractal limit. 

In the sequel, we will be interested in the "inertial 
range" which exists if £(2, y) is very slow function of the 
dimensionless variable y = r/i over the interval from 
some 2/0 > 1 to some 51 > t/o > 1. In this case and if 
£(2,2/) < 2, we will be entitled to speak about the mul- 
tifractal behavior of the random process 6X(t) for some 
wide range of scales and define the effective exponent 
£(2) as the minimal value 



C(2)=mmC(2,2/) 
y 



(25) 



as we derive below as an optimal and efficient definition. 
Having obtained C(2), we can obtain a first estimation 
of the intermittency coefficient A 2 through the relation 



A 2 = 2 - C(2) , 



(26) 



which assumes that the parabolic dependence (|10|) 
holds. 



C. Method of determination of the effective 
exponent £(2) 

As we pointed out, the absence of exact multifractility 
for ip > does not prevent the process from exhibiting 
approximate multifractal scaling which, for all practi- 
cal purposes, can be undistinguishable from an exact 
one. Indeed, empirical data is always noisy and power 
law scaling are always sampled on a finite (often small) 
range of scales j^, H3, 13 • The experimentally relevant 
question is thus to define scaling from an operational 
view point consistent with what is done empirically. We 
now describe a natural and robust determination of the 
approximate scaling which will be seen to link intrinsi- 
cally the definition and determination of the exponent 
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£(2) with the existence of an integral time scale L de- 
fined from the range over which the approximate power 
law scaling holds. 

For this, let us consider a y-interval [j/i, y?\. We define 
£(2, j/i, y 2 ) over this interval as 

ry2 

C(2,y 1)2 / 2 )=Arg min / dy\\\nS 2 {y)-{A+C\n{y)\\ 2 . 

(27) 

In words, expression l|27|) determines £(2,2/1,2/2) as the 
exponent £ which minimizes in a OLS sense the distance 
between In S 2 (y) and a straight line in the In y- variable 
over the interval [2/1 , 2/2] ■ Let us now introduce the pa- 
rameter 77, which measures the precision with which an 
approximate power law scaling is qualified. Specifically, 
a power law with exponent £(2, 2/1,2/2) is qualified if 

min/ dy\\\nS 2 (y)-(A + CHy)\\ 2 < V . (28) 
Jyi 

Otherwise, the power law scaling is rejected. For a fixed 
77, we scan all possible values of 2/1 and 2/2 for which con- 
dition (|28|l is verified and we select the couple (2/1,2/2) 
such that the range 2/2 / 2/1 is maximum. We thus obtain 
an approximate power law scaling with apparent expo- 
nent £(2,2/1,2/2) within the confidence or noise level 77 
over the maximum range [2/1 , 2/2] • We believe that this 
procedure embodies in a precise way the general fitting 
procedure of experimental data. 

Let us study the limit 77 — > 0. In the case where 
S 2 {y) is not a pure power law (and thus In £2 (2/) is not 
a perfect straight line In y), the condition Ij28(l imposes 
that 2/2/2/1 1- Therefore, £(2,2/1,2/2 — * 2/1) determined 
by H27fl yields the local slope of the function In 6*2(2/) i n 
the variable In 7/, in other words 

Consider now the function 

f yi . y2 ^yi(yi) = ln £2(2/1) ~ £(2,2/1,2/2 -> 2/l) hi(2/i) ■ 

(30) 

This function f yij y 2 ^ yi (2/1) has at least one minimum in 
the interval y G [£,+00], which we call 2/0 ■ Close to its 
minimum, the function f yi ,y 2 — >y A (2/1) can be expanded 
up to second order to obtain 

fviM-^in = ^ + (l/2)a(ln(yi) - Hvo)) 2 , (31) 
where 

a = d 2 \nS 2 /d 2 \n{y)\ yo (32) 

is the second-order derivative of ln 5*2 with respect to 
lnj/, estimated at j/o- It is clear that, for small fi- 
nite values of 77, the largest range 2/2/2/1 is obtained for 



2/1 = 2/0 (we assume here that / is convex so that the 
local minimum is the global minimum; for the range of 
2/'s in the examples we have investigated, we have found 
the convexity condition to always hold). In addition, 
/j/!,3/2— >yi (2/1) does not change appreciably over a range 
of ln(2/i) proportional to 1/a 1 / 2 . Thus, the smallest a 
is, the largest is the range over which f yi ,y 2 — >y 1 (yi) will 
be almost constant and thus over which the exponent 
£(2) will well-defined and constant. This reasoning pro- 
vides the algorithm to measure the approximate expo- 
nent £(2) for a given data set, which we are going to 
use in a systematic way. It is given by (I29|) , where 2/1 is 
chosen equal to the argument 2/0 such that the second 
derivative (|32[1 is zero (or reaches the minimum posi- 
tive value over the whole range available when zero is 
not crossed). This is equivalent to searching for the ex- 
ponent i|29[l which takes the smallest possible positive 
value over the range of study. 



IV. SCALING OF THE SECOND-ORDER 
MOMENT S 2 

We study the process whose properties are con- 
trolled by the two key parameters tp defined in Q and 
a 2 defined in l[TTJl. We calculate the second-order struc- 
ture function S 2 {y) defined in l|16|) with y = t/£, where 
G{x) is obtained from expressions (|13|l and i|14|) . 

As a first example, we fix <p — 0.01 and scan a 2 = 
20; 30; 40. The corresponding structure functions S 2 (y) 
are plotted in figure 1 as a function of y = t/£ in log- 
log scales, together with the best power law fits shown 
as straight lines. A superficial examination suggests an 
excellent scaling behavior over at least four orders of 
magnitudes, with an exponent which clearly varies with 
a 2 from 1.66 to 1.34 when a 2 goes from 20 to 40. This 
property is novel compared with the MRW and previ- 
ous multifractal process with tp = and, as we are going 
to explore in some length, allows us to control the mul- 
tifractal properties continuously as a function of a 2 as 
well as tp. Fig. 2 examines more precisely the nature 
of the apparent power law behavior depicted in Fig. 1 
by plotting the local exponent £(2,2/) defined in 118|) as 
a function of y in log-scale in the range y s [1,10 7 ], 
for a fixed a 2 = 10 and varying values of tp — 0.002 to 
0.01. The first important message of this figure is that 
the exponent £(2,2/) is approximately constant over a 
large range of y- values, all the more so, the closer <p is 
to zero (this later property is of course not a surprise 
since <p —> recovers previously known multifractal pro- 
cesses) . Interestingly, this approximately constant value 
for £(2,2/) has a rather large dependence on tp itself, 
showing again that we can control the intermittency 
parameter by changing tp as well as a 2 . Another im- 
portant observation is the sharp variation of £(2, y) on 
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the left-hand-side of the range, suggesting a rather well- 
defined "viscous scale" t v , which we characterize as the 
boundary between the region of approximate constancy 
of C(2, y) around the definition 1|25[1 and the region of 
sharp increase of £(2,y) as y decreases below the min- 
imum (|25ll . In the examples shown in fig. 1, t v is in 
the range 10 2 -j- 10 3 . The closer ip is to 0, the smaller 
is the viscous scale t v and the better is the multifrac- 
tal scaling. In contrast, it is not obvious to identify the 
integral scale over the interval shown in Fig. 1 as the 
increase of £(2, y) towards 2 has not yet occurred ap- 
preciably even up to y = 10 7 , ensuring a rather nice 
approximate scaling with anomalous multifractal expo- 
nent C(2) < 2. Another way to express this observation 
is that the "inertial regime" over which the scaling of 
the second-order moment holds has not absolutely de- 
fined boundaries. Therefore, the transition to the in- 
tegral scale is smooth, a property also documented in 
hydrodynamic turbulence [5j. 

Fig. 3 is the same as Fig. 1 but for a large value of 
tp = 0.5 and with a 2 — 1 and 5. It shows that the scaling 
law for 5*2 still holds over approximately two decades in 
the horizontal y scales (and of course more in the ver- 
tical scale). A plot analogous to Fig. 2 for these values 
of tp and a 2 shows obviously larger deviations from a 
constant behavior (not shown). As we discussed in the 
previous section, the important message is that we have 
a rather large latitude in changing tp and a 2 to exhibit a 
reasonable (with respect to standard experimental pre- 
cision) multifractal scaling. 

Fig. 4 shows the other limit of a small value of 
tp = 0.001 for a large range of values of a 2 = 100 
to 500. While we expect indeed that the multifractal 
scaling should extend on large ranges as ip decreases to 
zero, the most remarkable fact is that the exponent £(2) 
can be continuously adjusted at will from 2 all the way 
to close to 1 by varying a 2 at fixed tp. Note also the 
large range of y-scale over which the apparent exponent 
C(2, y) remains approximately constant: for instance, 
for a 2 = 100, tp = 0.001, we measure a well-defined con- 
stant exponent £(2) = 1.81 over 8 orders of magnitudes. 
Even for these very large ranges of scales, the known 
existence of an integral scale and the transition to the 
normal value £(2) = 2 is not seen. 

Fig. 5 gives a synopsis of the dependence of the ef- 
fective exponent ((2) as a function of the two control 
parameter tp and a 2 . Actually, we show instead the "in- 
termittent coefficient" A 2 = 2 — £(2) first introduced in 
equation (|2l)|> . in analogy with the parabolic multifrac- 
tal spectrum <|10|) . The most important observation is 
that A 2 increases linearly with a 2 before saturating to 1 
asymptotically for large <r 2 's. The upperbound 1 for A 2 
results from the following property of the second-order 



moment £2 obtained from the definition (| 1 fc>|) : 
d 2 S 2 (y) 



dy 2 



2G(y) > 



(33) 



The fact that the second-order derivative of S^y) is 
positive means that 5*2(2/) increases faster than a linear 
function of y, hence 



C(2) > 1 



A < 1 



(34) 



V. HIGHER-ORDER MOMENTS, UNIVERSAL 
SCALING FUNCTION AND MULTIFRACTAL 
SPECTRA 

A. Definition and determination of the 
higher-order moments 

We have also investigated the higher-order moments 
5*3, 5*4 and £5 up to order 5 of the increments 5 T X(t) 
defined in JSJ). Moments of arbitrary orders can be ob- 
tained from formulas generalizing expression (|16[) for S2 , 
as follows 

S q (y) = q(q - 1) / (y - x)G q {x)dx , (35) 
Jo 

where, for q > 2, 

C-X C-X 

G q (x) = G(x) / dui . . . / dll q -2 
Jo Jo 



q-2 

Yl G(xi)G{u - Xi)G(xi - Xj) 

i=l 
j=i+l 



(36) 



The corresponding local exponents for the higher-order 
moments, are defined analogously to i|18fl and (|19fl as 



C(q,y) 



2-A q (y) ' 



where 



xG q {x)dx 



G q {x)dx 



(37) 



(38) 



Using these relations, we adopt the definitions general- 
izing H25f) for the effective exponents 



C(q) = minC(g,?/) 

y 



(39) 



associated with the effective power laws 

S q (y) = A q y^ , A q = S q (y m )y m ^ , (40) 
where y m is the value of y which makes £(<7, y) minimum. 
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B. Universal scaling ansatz 

Our numerical calculations of the higher-order mo- 
ments {q ^ 5) show that excellent scaling is observed 
for these moments over a wide range of the dimension- 
less variable y, similarly to the case of the second-order 
moment S2 presented in figures 1-3 [25]. This allows 
us to determine the dependence of the effective multi- 
fractal spectrum £(g) defined by (|39[1 with respect to q, 
as well as a 2 and tp. For this, it is convenient to use 
the parabolic spectrum (|10[1 as a proxy to extract an 
effective (a priori) g-dependent intermittent coefficient 
defined by 



A 2 (<?; cr 2 ,^) 



, g-C(g) 

'<?(<?-!) 



(41) 



We also make explicit in the notation A 2 (q; a 2 , tp) the 
dependence on the two parameters a 2 and (p. Similarly 
to the case of the second-order moment, expression (13511 
allows us to show that ((q) and X 2 (q; a 2 , tp) satisfy to 
the following inequalities 



C(q) > 1 



X 2 (q;a 2 ,p) <2/q , 



(42) 



which generalizes (13411 for arbitrary g's. 

Our detailed numerical calculations suggest the con- 
jecture that the dependence of X 2 (q; a 2 , ip) with respect 
to cr 2 ,ip can be factorized as follows: X 2 (q; a 2 , tp) = 
2A(aq)/q, where the factor a — a(a 2 ,p) depends only 
on the parameters a 2 and tp and not on q and the func- 
tion A(x) is a monotonously increasing. Moreover, the 
numerical calculations of the intermittent coefficients 
show that it is an excellent approximation to represent 
a(a 2 ,ip) as a linear function of a 2 , i.e., a = b(ip)cr 2 . This 
provides the following universal scaling law for the gen- 
eralized intermittency coefficients 



X 2 (q;<J 2 ,v) = -A(ba 2 q) 

q 



(43) 



The factor b, which is independent of q and of a 2 accord- 
ing to the scaling ansatz (|43|) . can be determined from 
the intermittency coefficient X 2 (q = 2; a 2 , tp) = A(2ba 2 ) 
for q = 2, previously reported. Since X 2 (q — 2;a 2 ,tp) is 
a linear function of a 2 for not too large <r 2 's as shown in 
Fig. 5, this implies that the function A(x) ~ x is linear 
for small x J s. Defining b such that, for x — > 0, we have 
A(x) ~ x, we obtain 



%0 



A 2 (g = 2;a 2 ,y>) 
2a 2 



(A 2 « 1) . (44) 



Taking for instance a 2 = 20 for which all intermittency 
coefficients remain small, we find that b(ip) is very well 
represented by the following power law dependence 



Fig. 6 plots the reconstructed functions A(x) = 
(q / 2) X 2 (q; a 2 , ip) for different orders q. The excellent 
collapse provides a first validation of the scaling ansatz 
I®. 



C. Multifractal spectra 

The general scaling ansatz <|43|) allows us to express 
the multifractal spectrum under the following form 



C(q) = q+(l~q)A(ba 2 q) 



(46) 



Moreover if, for some particular multifractal phe- 
nomenon, the intermittency coefficient A 2 = A 2 (2) is 
small compared to 1, then one can rewrite expression 
H46|) in the more universal form 



C(?) 



(1 - q)A(X 2 q/2) 



(47) 



oupf 



a = 0.58, /? = 0.92 



(45) 



whose dependence on a 2 and ip is completely embed- 
ded in that of A 2 . This implies that the knowledge of 
the intermittency coefficient A 2 together with the scaling 
function A(x) allows one to determine the multifractal 
spectrum for arbitrary orders even for q 1 . 

Our prediction (|47|l can be checked by direct numer- 
ical calculations of the moments up to order 5 that we 
have performed. Fig. 7 plots the multifractal spectra 
((q) for q = to 5 obtained by two methods: (i) the 
circles are the direct numerical integration of (|35l3ti|) of 
the stochastic process for tp — 0.004 and different val- 
ues of a 2 ; (ii) the continuous lines are obtained by using 
Ij46(l with the scaling function A(x) constructed as in 
Fig. 6 for q = 2 and ip — 0.001 and applied to the case 
tp = 0.004. The agreement is good, even for the large 
values q = 5. This illustrates that different combina- 
tions of ip and a 2 with fixed ba 2 give the same value of 
the intermittency coefficient A 2 (for A 2 not too large) ac- 
cording to (|43l) and thus the same multifractal spectrum 

C(?)- 

Note that £(g) becomes non-concave for large tr 2 's at 
large g's, a property which is excluded for exact multi- 
fractal scaling by the Holder inequality applied to the 
moments {[6 T X(t)] q ). The absence of concavity in a 
certain range of parameters reflect the fact that mul- 
tifractality is not an asymptotic property observed at 
small or large values of the dimensionless variable y, but 
only in some intermediate scaling range. Non-concavity 
also prevents obtaining the exact multifractal spectrum 
of dimensions f(a) of singularities a, but yields only a 
concave envelop of it j2(| . 
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VI. NUMERICAL AND MATHEMATICAL 
INSIGHTS ON THE ORIGIN OF THE 
EFFECTIVE MULTIFRACTALITY 

As we showed, these multifractal properties are only 
effective properties or approximations, as exact multi- 
fractality only holds when C(r) defined in l|14f> is pro- 
portional to the logarithm of r (up to an integral scale). 
The approximate multifractality discussed here can be 
tracked back to the approximate logarithmic depen- 
dence of C(t), as we know make clear, numerically and 
mathematically. 

Let us consider the second-order moment S^y). The 
scaling l(TT}l is a precise description of 5*2 (y) if G(y) in 
(I13|) is close to a power law, i.e., a 2 d(y) is close to A 2 Iny. 
To test if this is the case for different values of a 2 and 
ip, we construct 



S(y) = a 2 d(y)/\ 2 



(48) 



which should be close to In y to justify our results above. 
Recall that, for rather small a 2 (actually, we just need 
that a 2 < 100), we have the relation lO with igSJ, 
which together with l|48|) yields 



d(y) _ 1 - C(y) 

K<p) Hfp) 



(49) 



which has to be almost equal to In y for any a 2 and ip to 
justify our results. This is verified in Fig. 8, which plots 
S(y) given by l|49[) . As ip departs more and more from 0, 
increasing deviations from lny occur, which reduce the 
range over which scaling and multifractality can be ob- 
served. An argument similar to that presented in Fig. 8 
was proposed in 0, ^| for the multi-scaling of the con- 
ditional response function (Omori law). 

This argument can be made mathematically precise 
as follows. Let us consider the function 



K(y,<p) = / h(x,ip)h(x + y,cp)dx 



where we re-define 



h(x,ip) = 



such that 



<j* = K{0,<p) = — 



(50) 



(51) 



(52) 



Determining how well the covariance function C(t), de- 
fined by H14|). is approximated by a logarithm of t 
amounts to study how well K(y 1 ip) is approximated by 
a logarithmic function of y. 



The study of K(y, ip) is based on its explicit analytical 
expression 



K(y,(p) 
-2 



l-2ip 



/tt T(<p) sec(ir<p) f 2\ 2ip 

2IV+T72) \y) (53) 
-F (\,\+<p,\-tp,\ + y) , 



where F(a, b, c, z) is a hypergeometric function. 

Let us introduce the analog of the function d(y) (see 
O) defined by 

D(y,ip)=a 2 -K(y,ip) = ^--K{y,ip) . (54) 

zip 

Using H53|) . we obtain the following asymptotic expres- 
sion for D(y, ip): 



D(y,f) 



l F(y,ip)-F(y,0) 
2 ip 



In 



VTT^-l 



valid for — > 0, where 



F(y,<p) = p(<p)[- 

\y 



2(p 



and 



p(v>) = ^(99) 



m/2+ip) 



. sec(irip) 



(55) 
(56) 

(57) 



We verify that the asymptotic expression (|55|l gives 
an excellent approximation to the exact expression for 
D(y, ip) obtained by using (|53|l for 1 < y < 10 15 for 

To reveal the dependence of D(y, ip) with respect to 
y for small tp, we may interpret the first fraction in the 
r.h.s. of the relation (|55|l as a discrete approximation 
of the derivative F'(y,(p) of the function F(y,ip>) with 
respect to tp. This leads to the following approximate 
relation 



D{y,v)' 
where 

F\y,v) 

with 

A(ip) 



H4) 



In 



/yTT^+1 



-2 In 



A((p)y 



p(v) 



exp 



(58) 



(59) 



(60) 



Since In [(yT+y + l) / (yT+y - 1) ] = 2/^ + 
0(1/ y) for y ^> 1, we can neglect this logarithmic term 
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in (|58|l for p < 1/2 and obtain the following logarithmic- 
power approximation 



D(y,<p)~ hi A(p/2)y^ 



2 \ * 



(61) 



Figure 9 plots the exact function D(y, ip) and its ap- 
proximation 1|61|) . showing the approximate linear de- 
pendence of D(y,p) as a function of hiy for a large 
range of y for ip > 0. 

For ip — > 0, p(p/2) — > 1 and we recover D(y,p — > 
0) — > lny, as expected. For y> > 0, expression (ffTTft 
makes explicit the deviation from a pure logarithmic 
dependence of -D(y, </?) in the form of the power factor 
{2/y) v . The predicted deviation from a pure logarithm 
gives us the possibility to estimate the integral scale 
beyond which the effective multifractality breaks down. 
We thus define the integral scale L e as the solution of 
the equation 



2 X ¥> 

V, 



(62) 



where e < 1 measures the deviation of D(y, tp) from lny. 
This gives L £ = 2/e 1 /^. Taking arbitrarily e = 1/2, we 
obtain L 1/2 ~ 2000 for tp = 0.1, L 1/2 - 1.3 x 10 30 for 
p = 0.01 and the astronomical value Li/ 2 ~ 10 301 for 
tp = 0.001. These values explain why the saturation of 
the effective multifractal scaling predicted for large y is 
not observed for small p's. 



VII. CONCLUDING REMARKS 

We have confirmed on the multifractal spectrum 
the proposition previousl y in troduced in the context of a 
model of earthquake [la, that processes constructed 



as exponentials of long-memory processes should exhibit 
multifractal properties over a significant range of the 
parameters. While the initial argument concerned the 
time decay of the conditional expectation of the response 
function (called the Omori law for the decay of the num- 
ber of aftershocks in the context of seismology) as a 
generalization to the prediction and observation for the 
multifractal random walk discussed in Ref. [l5j , we have 
extended the analysis to show that multifractality is a 
robust property of this class of processes defined as ex- 
ponentials of long memory processes. 

Notwithstanding the fact that the multifractal prop- 
erties discusses here are only effective properties or ap- 
proximations, we claim that there is probably no way to 
tell the difference between an exact multifractal random 
walk or an exact multifractal scaling from our approxi- 
mate scaling laws and approximate £(q) functions, given 
the noise and the limited ranges usually seen in experi- 
mental and numerical studies. We have illustrated this 
claim by showing that a given multifractal spectrum 
can be obtained for several sets of parameters of our 
model. In a sequel, we will discuss what observables 
should be used to possibly distinguish between different 
implementations, so as to extract useful physical inter- 
pretation of the parameters A 2 , p and a 2 . 
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Fig. 1: Log-log-plot of the second-order mo- 
ment 5*2 (y) 1161 1 and its power approximation 
(fT7|l for ip = 0.01 and a 2 = 20; 30; 40 (top to bot- 
tom). The corresponding exponents are equal to 
C(2) = 1.66; 1.49; 1.34. 




Fig. 2: Local exponent C,{2,y) as a func- 
tion of y = t/£ for u — 10 and ip = 
0.002; 0.004; 0.006; 0.008; 0.01 (top to bottom). 
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Fig. 3: Log-log plot of second-order moment 
S2(y) calculated using 1161 and its power ap- 
proximation ifY7\ for ip = 0.5 and a 2 — 1; 5 (top 
to bottom). The corresponding effective expo- 
nents are respectively £(2) = 1.82; 1.26. 
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Fig. 4: Dependence of the local exponent 
f(2, y) as a function of y = t/1 in log-scale, for 
ip = 0.001 and a 2 = 100; 200; 300; 400; 500 (top 
to bottom). 
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a 2 

Fig. 5: Dependance of the intermittency co- 
efficient A 2 = 2 — f(2) as a function of a 2 for 
different values of ip = 0.01 -j- 0.04 (bottom to 
top). 







0.5 

x 



1 



1.5 



Fig. 6: Plot of the universal scaling function 
A(a;), obtained from relations (1411 and 143 I I and 
the numerical calculation of the dependence of 
the effective multifractal exponents £(q) as a 
function of a 2 , for <p — 0.001 and q — 2; 3; 4. 
The slight discrepancies between the curves in 
the neighborhood of x = 1 can be attributed 
to some systematic errors of numerical calcula- 
tions. 
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Fig. 7: Universal multifractal spectra ((q) for 
tp = 0.004 and a 2 = 10; 20; 30; 40; 50; 60 (top 
to bottom). We compare two different methods 
for estimating C(q): (i) the circles are the direct 
numerical integration of 1351361 ; (ii) the contin- 
uous lines are obtained by using l|46|l with the 
scaling function A(x) constructed as in Fig. 6 
for q = 2. 
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Fig. 8: Dependence of S(y) given by l|49|l 
as a function of y in logarithmic scale for 
ip = 0.01,0.02,0.05,0.1 (top to bottom). The 
straight dashed line corresponds to the logarith- 
mic dependence \ny — 0.8. Not surprisingly, the 
closer ip is to 0, the larger is the range of lny 
over which the approximation 5(y) ~ \ny holds. 



Fig. 9: Plots of the exact function D(y, ip) and 
its approximation l|61|l for ip = 0.001,0.01,0.02 
(top to bottom) 
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